On the instability of Reissner-Nordstrom black holes in de Sitter backgrounds 
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Recent numerical investigations have uncovered a surprising result: Reissner-Nordstrom-de Sitter 
black holes are unstable for spacetime dimensions larger than 6. Here we prove the existence of 
such instability analytically, and we compute the timescale in the near-extremal limit. We find very 
good agreement with the previous numerical results. Our results may me helpful in shedding some 
light on the nature of the instability. 
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I. INTRODUCTION 

In physics, stability of a given configuration (solution 
of some set of equations), is a useful criterium for rele- 
vance of that solution. Unstable configurations are likely 
not to be realizable in practice, and represent an interme- 
diate stage in the evolution of the system. Nevertheless, 
the instability itself is of great interest, since an under- 
standing of the mechanism behind it may help one to 
better grasp the physics involved. In particular, it is of 
interest to be able to predict which other systems display 
similar instabilities, or even have a deeper understanding 
of the physics behind the instability (why is the system 
unstable? is there some fundamental principle behind 
the instability?). 

In General Relativity, the Kerr family exhausts the 
black hole solutions to the electro- vac Einstein equations. 
Kerr black holes are stable, and can therefore describe 
astrophysical objects. However, there are many instances 
of instabilities afflicting objects with an event horizon, 
such as the Gregory-Laflamme the ultra-spinning ^] 
or superradiant instabilities 3 and other instabilities of 
higher-dimensional black holes in alternative theories [1, 
Q(for a review see Ref. [6|). 

Konoplya and Zhidenko (hereafter KZ) recently stud- 
ied small perturbations in the vicinity of a charged black 
hole in de Sitter background, a Reissner-Nordstrom de 
Sitter black hole (RNdS) 0. Their (numerical) results 
show that when the spacetime dimensionality D > 6, the 
spacetime is unstable, provided the charge is larger than 
a given threshold, determined by KZ for each D. Because 
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the results are so surprising (the mechanism behind it is 
not yet understood), we set out to to investigate this in- 
stability and hopefully understand it better. Our results 
can be summarized as follows: (i) we can prove analyti- 
cally the existence of unstable modes for charge Q higher 
than a certain threshold, (ii) in the near-extremal regime, 
we are able to find an explicit solution for the unstable 
modes, determining the instability timescale analytically. 
We hope that our incursion in this topic helps to better 
understand the physics at work. 



II. EQUATIONS 

This work focuses on the higher dimensional RNdS ge- 
ometry, described by the line element 



-fdt^ + f-^dr^+r''dn 



where dfl^ is the line element of the n sphere and 



/ = 1-Ar2 



2M 



(1) 



(2) 



the background electric field is Eq = q/r", with q the 
electric charge. The quantities M and Q are related to 
the physical mass M and charge q of the black hole [1], 
and A to the cosmological constant. The spacetime di- 
mensionality is D — n 2. 

The above geometry possesses three horizons: the 
black-hole Cauchy horizon at r = , the black hole event 
horizon is at r = rt and the cosmological horizon is at 
r = Tc, where rc > r;, > Tq, the only real, positive zeroes 
of /. For convenience, we set rh = 1, i.e., we measure all 
quantities in terms of the event horizon rf,. We thus get 
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2M^1 + Q^ -X, 



(3) 
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Furthermore, we can also write 



A 



frl) 



r, 



n+2 



(4) 



For a fixed Tc and spacetime dimension D, the existence 
of a regular event horizon imposes that the charge Q must 
be smaller than a certain value Qext- With our units this 
maximum charge is 



Qoxt ~ 



{-2rc + {n+ l)r? - (n - l)r?+2) 
-fc (?'c(?i + 1) - 2nr^' + (n - l)r?"+^) 



(5) 



Gravitational perturbations of this spacetime couple to 
the electromagnetic field, and were completely character- 
ized by Kodama and Ishibashi 8] . They can be reduced 
to a set of two second order ordinary differential equa- 
tions of the form, 



where the tortoise coordinate and the potentials Vs± 
are defined through 



We have 



fU± 



2 

n{n + 1) 



(1 + mS)x , 



(7) 



(8) 
(9) 



Vs±) $± - , 



(6) and the quantities U± are given by 



U+ = [-4n^{n + 2){n + l)^5^x^ - 4:8n^{n + l){n-2)Sx 

-16{n - 2){n - 4)] y - S^n^{3n - 2)(n + 1)'*(1 + mS)x'^ 

+4:5^n^{n + i f {{n + l){3n - 2)m5 + in'^ + n ~ 2} x^ 

+4:5{n + 1) {{n - 2){n - 4)(n + l)(m + n^K)S - 7n^ + Ir? - 14n -f 8} 

-h{l6(n+l) (-4m + 3n2(n-2)X)(5- 16(3n - 2)(n - 2)} a; 

+64to + 16n(n 2)if, (10) 
[/_ = [-4n3(n + 2)(n-h 1)^(1 + m(5)2a;2 + 48n2(n+l)(n-2)TO(H-TO(5)a; 
-16(n - 2){n - 4)™^] y _ ^3^3^ _ 2)(n + \f8{\ + m8fx^ 
-\n^{n + 1)2(1 + m6f {[n + l)(37i - 2)m5 - n^) 
+4(n+ 1)(1 + m(5) {m{n - 2)(n - 4)(n + l)(m -|- n2i4:)(5 
-f4n(2n2 - 3n + 4)m + ^^(n - 2)(n - 4)(n + 1)^} a;^ 
-16to {(n l)m (-4m + Zr?(n - 2)^) (5 
-h3n(n - 4)m + 3n2(n l)(n - 2)is:} a; 

-h64m^ -f 16n(n + 2)m2i4:. (11) 



The variables x, y and parameters fi, m are defined 
through 



_ 2M 
X = , y = Xr 



(n + l)2 



(12) 

■m = k^-nK, (13) 



and the quantity 5 is imphcitly given by /x = (l + 2m(5)M. 
Note that the following relations holds — (71 + 
lfM^5{l + m5). 

Note also that for the spacetime considered in this pa- 
per K = \, which means that the eigenvalues fc^ are given 



by fc^ = l{l + n — 1), where I is the angular quantum 
number, that gives the multipolarity of the field. The 
behavior of the potentials varies considerably over the 
range of parameters. In Fig. [1] we show y_ for = 8, 
Tc = 1/0.95, ? = 2 and three different values of the charge, 
Q = 0.2,0.35,0.44. 



III. A CRITERIUM FOR INSTABILITY 

A sufficient (but not necessary) condition for the exis- 
tence of an unstable mode has been proven by Buell and 
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FIG. 1: Behavior of V- for different parameters, for D — 8. 
Here we fix the event horizon at r;, = 1, and the cosmological 
horizon at rc = 1/0.95. We consider I — 2 modes and three 
different charges, Q = 0.2, 0.35, 0.44. 



Shadwick [9|] and is the following, 

V 



f 



dr <0. 



(14) 



The instability region is depicted in figure [2] for several 




rb/rc 



FIG. 2: The parametric region of instability in Q/Qcxt 
coordinates, according to criterim (|14[) . for / — 2. 
bottom, D = 7,8,9, 10, 11. 



Top to 



spacetime-dimension D, which can be compared with the 
numerical results by KZ, their figure 4. It is apparent 
that condition (jl4p very accurately describes the numer- 
ical results for rh/vc ~ 1, a regime we explore below in 
Section IIVI As one moves away from extremality cri- 
terium p4p is just too restrictive. An improved analysis 
and refined criterium would be necessary to describe the 
whole range of the numerical results. Nevertheless, figure 



[2] is very clear: higher-dimensional {D > 6) RNdS black 
holes are unstable for a wide range of parameters. 



IV. AN EXACT SOLUTION IN THE NEAR 
EXTREMAL RNDS BLACK HOLE 

Let us now specialize to the near extremal RNdS black 
hole, which we define as the spacetime for which the cos- 
mological horizon rc is very close (in the r coordinate) 
to the black hole horizon rb, i.e. ^-^^f^ ^ 1- The wave 
equation in this spacetime can be solved exactly, in terms 
of hypergeometric functions . The key point is that 
the physical region of interest (where the boundary con- 
ditions are imposed), lies between rb and rc- Thus, 



(r - rb){rc - r) 



n 



(15) 



where we have introduced the surface gravity Kb associ- 
ated with the event horizon at r = rb, as defined by the 
relation Kb = ^df /drr^rt- For near-extremal black holes, 
it is approximately 



Kb 



(rc - rb){n - 1) 



2rl 



(16) 



In this limit, one can invert the relation r*(r) of ([7]) to 
get 



1 



(17) 



Substituting this on the expression (|15p for / we find 
(rc - rb)Kb 



f 



2cosh(Khr*) 



(18) 



As such, and taking into account the functional form of 
the potentials for wave propagation, we see that for the 
near extremal RNdS black hole the wave equation ^ is 
of the form 



d'^^{uj, r) 



dr'i 



cosh (Kbr^,) 



$(u;,r)=0, (19) 



with 



{rc ~ rb)Kb Vs±{rb) 
2 / 



(20) 



The potential in ([TO)) is the well known Poshl- Teller po- 
tential [ni- The solutions to pJ)) were studied and they 
are of the hypergeometric type, (for details see Refs. 
[l^fTsj). It should be solved under appropriate boundary 
conditions: 



— 00 

oo. 



(21) 
(22) 
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These boundary conditions impose a non-trivial condi- 
tion on uj [13, [13] , and those that satisfy both simultane- 
ously are called quasinormal freq uencies. For the Poshl- 
Teller potential one can show jl2l . ITsI that they are given 

by 




J =0,1, 



(23) 

We conclude therefore that an instability is present 



TABLE I: The threshold of instability for near-extremal 
RNdS black holes (i.e., black holes for which the cosmological 
and event horizon almost coincide) for / — 2 modes. We show 
the prediction from the exact, analytic expression obtained 
in the near extremal limit (|24p . which we label Q/Q^^t and 
the one from criterium (|14p which we label as Q/QY^t- Both 
these results are compared to the numerical results by KZ. 

D 



10 



11 



Q/Qcxt 0.913 0.774 0.683 0.617 0.567 V27d 



Q/Q. 

Q/Qi 



V 
cxt 




0.913 0.775 0.684 0.618 0.568 



0.94 0.78 0.68 0.61 0.55 — 



whenever Vq is negative. The threshold of stability in 
the near-extremal regime is therefore given by 



Vs±{rb) 
f 



0, 



(24) 



The expression for Vs±{rb) / f is lengthy, and we won't 
present it here. The values of the charge Q/Qext that sat- 
isfy the condition above are given in Table I (for 1 — 2), 
and compared to the prediction from the analysis in Sec- 
tion lllli criterium (fT^ . The agreement is excellent. Fur- 



thermore, we compare these predictions against the nu- 
merical results by KZ, extrapolated to the extremal limit 
(p = 1 in KZ notation). The agreement is remarkable. 



V. CONCLUSIONS 

We have shown analytically that charged black holes 
in de Sitter backgrounds are unstable for a wide range of 
charge and mass of the black hole, confirming previous 
numerical studies The stability properties of the ex- 
tremal D = 6 black hole remain unknown. Our methods 
and results and inconclusive at this precise point, further 
dedicated investigations would be necessary. 

Our analytical result in the near-extremal regime could 
be used to investigate further the nature of this instabil- 
ity, something we have not attempted to do here. A 
possible refinement concerns the large- Z? limit of the in- 
stability, where it could be possible to find an analytical 
expression throughout all range of parameters. We have 
in mind results and techniques similar to those of Kol and 
Sorkin It would also be interesting to investigate 

the stability properties, using this or other techniques, of 
near-extremal Kerr-dS black holes, which have recentl;v 
been conjectured to have an holographic description [1^ 
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